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The minimum bias multiplicity distribution and the transverse momentum and pseudorapidity 
distributions for central collisions have been measured for negative hadrons (h'^) in Au+Au in- 
teractions at ^Snn ~ 130 GeV. The multiplicity density at midrapidity for the 5% most central 
interactions is dNf^- / dri\^=o = 280 ± l(stat) ± 20(syst), an increase per participant of 38% relative 
to pp collisions at the same energy. The mean transverse momentum is 0.508 ± 0.012 GeV/c and is 
larger than in central Pb+Pb collisions at lower energies. The scaling of the yield per participant 
is a strong function of p_L. The pseudorapidity distribution is almost constant within \rj\ < 1. 



The collision of high energy heavy ions is a promising 
laboratory for the study of nuclear matter at high energy 
density and the possible creation and decay of the Quark 
Gluon Plasma . A new era in the study of high energy 
nuclear collisions began in the year 2000 with the first 
operation of the Relativistic Heavy Ion Collider (RHIC) 
at Brookhaven National Laboratory. 

The multiplicity and inclusive single particle transverse 
momentum {p±) distributions of hadrons are important 
tools for understanding the evolutionary path of the sys- 
tem created in the collision and help to determine the 
characteristics of the early, hot and dense phase. In this 
Letter we present the minimum-bias multiplicity distri- 
bution and the pseudorapidity (rj) and p±^ distributions 
for the 5% most central collisions for negative hadrons 
{h^) in Au+Au interactions at a center-of-mass energy 
of V^NN — 130 GeV per nucleon pair, measured by the 
STAR detector at RHIC. The results are compared with 
reference data from pp collisions at a similar energy and 
collisions of heavy nuclei at a lower energy. 

The main tracking detector for STAR is a large Time 
Projection Chamber (TPC), which measures charged 
particles in the pseudorapidity range \r]\ < 1.8 with com- 
plete azimuthal acceptance. It is placed inside a uni- 
form solenoidal magnetic field of strength 0.25 T. The 
trigger detectors are an array of scintillator slats (CTB) 
arranged in a barrel surrounding the TPC, and two 



hadronic calorimeters (ZDCs) at ±18 m from the detec- 
tor center and at zero degrees relative to the beam axis. 
The ZDCs intercept spectator neutrons from the collision 
and provide a measure of the collision centrality. Further 
details on the apparatus can be found in 

During the Summer 2000 run, RHIC delivered colli- 
sions between Au nuclei at y/s^^ = 130 GeV. The data 
presented here are from a minimum-bias sample, trig- 
gered by a coincidence of signals above threshold in both 
ZDCs with the RHIC beam crossing. The ZDC thresh- 
old was set to ensure efficient detection of single spectator 
neutrons. The efficiency of the ZDC coincidence trigger 
for central events was measured using a high-threshold 
CTB trigger. The trigger efficiency was found to be above 
99% for the entire range of multiplicities reported in this 
Letter. 

The offline reconstruction found a primary vertex for 
each event by propagating the measured tracks through 
the field towards the beamline. The vertex resolution for 
high multiplicity events is approximately 150 /im, both 
perpendicular and parallel to the beam axis. Events used 
in the analysis have a vertex within ±95 cm of the center 
of the TPC along the beam axis. The vertex finding ef- 
ficiency is 100% for events with more than 50 primary 
tracks in the TPC acceptance, decreasing to 60% for 
those with fewer than 5 primary tracks. 60,000 minimum- 
bias Au±Au events were used for this analysis. 
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Particle production was studied through the yield of 
primary negative hadrons, comprising mostly tt" with 
an admixture of and p. The h~ distribution in- 
cludes the products of strong and electromagnetic de- 
cays. Negatively charged hadrons were studied in order 
to exclude effects due to participant nucleons. Charged 
particle tracks reconstructed in the TPC were accepted 
for this analysis if they fulfilled requirements on number 
of points on the track and on pointing accuracy to the 
event vertex. The measured raw distributions were cor- 
rected for acceptance, track reconstruction efficiency, and 
contamination due to interactions in material, misidenti- 
fied non-hadrons, and the products of weak decays. The 
reconstruction efficiency was determined by embedding 
simulated tracks into real events at the raw data level, 
reconstructing the full events, and comparing the simu- 
lated input to the reconstructed output. This technique 
requires a precise simulation of isolated single tracks, 
achieved by a detailed simulation of the STAR apparatus 
based on GEANT [| and a microscopic simulation of the 
TPC response. The multiplicity of the embedded tracks 
was limited to 5% of the multiplicity of the real event 
in the same phase space as the simulated data, thereby 
perturbing the real event at a level below the statistical 
fluctuations within the event sample. 

The acceptance is on average 95% for tracks within 
the fiducial volume having p±^ > 300 MeV/c. The track- 
ing efficiency ranges between 70%-95%, depending on 
p± and the multiplicity of the event. For tracks with 
p± > 200 MeV/c the efficiency is above 85%. Accepted 
tracks for this analysis have 0.1 < p± < 2 GeV/c and 
\r,\ < 1.0. 

Instrumental backgrounds due to photon conversions 
and secondary interactions with detector material were 
estimated using the detector response simulations men- 
tioned above, together with events generated by the HI- 
JING model [Q. The simulations were calibrated us- 
ing data in regions where background processes could 
be directly identified. The measured yield also contains 
contributions from the products of weak decays, primar- 
ily ATg, that were incorrectly reconstructed as primary 
tracks. The background fraction of the raw signal is ap- 
proximately 20% at p±= 100 MeV/c, decreasing with in- 
creasing p± . The average fraction of background tracks 
in the uncorrected sample is 7%. All corrections were 
calculated as a function of the uncorrected event mul- 
tiplicity. The systematic uncertainty due to the correc- 
tions was estimated by studying the variation in the final 
spectra due to both a large variation in the track quality 
cuts with corresponding recalculation of the correction 
factors, and a small variation in the track quality cuts 
with no adjustment of the correction factors. 

Figure ^ shows the corrected, normalized multiplicity 
distribution within \ri\ < 0.5 and p± > 100 MeV/c for 
minimum bias Au-t-Au collisions. The data were normal- 
ized assuming a total hadronic inelastic cross section of 



7.2 barn for Au-|-Au colfisions at ^/s^= 130 GeV, de- 
rived from Glauber model calculations Q . The multiplic- 
ity bin below Ni^- = 5 is not shown, due to large system- 
atic uncertainties in the vertex reconstruction efficiency 
and large background contamination. Its relative contri- 
bution to the total cross section was estimated to be 21% 
by normalizing the HIJING multiplicity distribution to 
the measured data in the region 5 < N/^- < 25. This 
procedure relies on the assumption that very peripheral 
interactions are well described by the superposition of a 
few nucleon-nucleon collisions in the geometry of a nu- 
clear collision, and can therefore be accurately modelled 
by HIJING. The systematic error on the vertical scale 
is estimated to be 10% and is dominated by uncertain- 
ties in the total hadronic cross section and the relative 
contribution of the first bin. The systematic error on the 
horizontal scale is 6% for the entire range of multiplicity 
and is depicted by horizontal error bars on a few data 
points only. 

The shape of the multiplicity distribution is domi- 
nated over much of the Nf^- range by the nucleus-nucleus 
collision geometry, consistent with findings at lower ener- 
gies. However, the shape of the tail region at large Nf^- is 
determined by fluctuations and acceptance. These over- 
all features are also observed in the HIJING calculation, 
shown as histogram in Fig. ^ The distribution for the 
5% most central collisions (360 mb), defined using ZDC 
coincidence, is shown as the shaded area in Fig. |l|. 

Figure |^, upper panel, shows the transverse momen- 
tum distribution of negatively charged hadrons for cen- 
tral Au-|-Au collisions at midrapidity {\r]\ < 0.1) within 
0.1 < p± < 2 GeV/c. Statistical errors are smaller than 
the symbols. The correlated systematic error is esti- 
mated to be below 6%. The data are fit in the range 
0.2 < p± < 2 GeV/c by a QCD inspired power-law 
function of the form d?Ni^- / dp\drj ~ A{1 +P_l/po) " 
where n, and po are free parameters. The upper 
panel of Fig. || also shows the pj^-distributions of neg- 
atively charged hadrons for central Pb+Pb collisions at 
y/s^— 17 GeV from NA49 and for minimum-bias 
pp coflisions at 200 GeV from UAl 0, fitted with 

the same function. The NA49 distribution, which was 
reported in units of pion rapidity, was transformed to 
units of pseudorapidity. The UAl invariant cross section 
Ed?'a/d?'p reported in Ref. was scaled by 27r/cri,iei, 
where (Tinci — 42 mb The power law fits all three 
datasets well. The mean can be derived from the fit 
parameters as (p^) = 2po/{n — 3). The fit to the STAR 
data gives pa = 3.0 ± 0.3 GeV/c, n = 14.8 ± 1.2, and 
(p^) = 0.508 ± 0.012 GeV/c. The strong correlation of 
fit parameters po and n must be taken into account when 
calculating the error on The (p^) from STAR is 

larger than that from both central collisions of heavy nu- 
clei at much lower energy ((p±)na49 — 0-429 GeV/c ||]) 
and nucleon-nucleon collisions at a comparable energy 
((p±)uAi = 0-392 ± 0.003 GeV/c). 
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The PHENIX collaboration has reported that 
{E±) / (Nch) for central collisions of heavy nuclei is con- 
stant between SPS and RHIC energies, with a value of 

0.8 GeV Under the assumption that the particle 
composition does not change significantly between the 
SPS and RHIC, this finding is in apparent disagreement 
with the above observation that {p±) increases by 18%. 
It should be noted, however, that the systematic error on 
the SPS measurement of {E±_)/{Nch) is 20% |ll), so that 
these results are consistent within errors. 

Figure ^ lower panel, shows the ratio of the STAR 
and UAl pj_-distributions. Since the UAl distribution 
is measured at y/s= 200 GeV, da/dp± is scaled by two 
factors for quantitative comparison to the STAR data at 
130 GeV: (i) i?(130/200), the pj^-dependent ratio of in- 
variant cross sections for charged particle production in 
pp colHsions at \/s— 130 and 200 GeV, and (ii) Taa= 
26 ± 2 mb^^, the nuclear overlap integral ||l^ for the 
5% most central Au-|-Au collisions. R varies from 0.92 
at p± = 0.2 GeV/c to 0.70 at pi_ = 2.0 GeV/c, and 
was derived using scaling laws for {pj_) and dNch/drj as a 
function of ^/s |0,|l^ together with the extrapolation to 
130 GeV of power-law parameterizations at = 200- 
900 GeV 0. The shaded boxes show the total error of 
the ratio, which is the linear sum of the errors of the 
measured data, depicted by the error bars, and the sys- 
tematic error due to uncertainties in the scaling with Taa 
and R. 

There are two simple predictions for the scaled ratio. 
In lower energy hadronic and nuclear collisions, the total 
pion yield due to soft (low p±) processes scales as the 
number of participants ( "wounded" nucleons) in the col- 
lision (e.g. The scaled ratio in this case is 0.164, 
assuming 172 participant pairs p5[ | and a mean num- 
ber of binary collisions of 1050 (= (JinciTAA, where for 
= 130 GeV (Tinci ^ 40.5 mb g) for the 5% most cen- 
tral Au-(-Au events. In contrast, if hadron production is 
due to hard (high p\^ processes and there are no nuclear- 
specific effects (see below), the hadron yield will scale as 
the number of binary nucleon-nucleon interactions in the 
nuclear collision and the value of the ratio is unity. There 
are important nuclear effects which may alter the scaling 
as a function of p±^ from these simple predictions, includ- 
ing initial state multiple scattering |16|, shadowing 
jet quenching jl^, and radial flow |1£]. The scaled ra- 
tio exhibits a strong pj^ dependence, rising monotonically 
with increasing p^ from Wounded Nucleon scaling at low 
Pj^ but not reaching Binary Collision scaling at the high- 
est p± reported. This is consistent with the presence of 
radial flow, as well as the onset of hard scattering contri- 
butions and initial state multiple scattering with rising 

Figure |^ shows the normalized pseudorapidity distri- 
bution of /i^ for the 5% most central collisions within 
\-q\ < 1.0, both for p±_ > 100 MeV/c and for all p^. 
The latter was obtained by fitting a power-law function 



in the range 0.1 < p± < 2 GeV/c and extrapolating to 
p± = in order to estimate the content of the first bin. 
The error bars indicate the uncorrelated systematic er- 
rors. The statistical errors are negligible. The correlated 
systematic error applied to the overall normalization is 
estimated to be below 6% for p± > 100 MeV/c and 7% 
for all p± . 

The 7] distribution is almost constant within \r]\ < 1, 
exhibiting a small rise at larger r]. This shape is ex- 
pected from a boost invariant source (i.e., constant in 
rapidity), taking into account the transformation from 
y to rj. It should be noted that in Pb-|-Pb collisions 
at ^/s^= 17 GeV the pseudo-rapidity distribution of 
charged hadrons |Q and the rapidity distribution of neg- 
ative hadrons (assuming the pion mass) were found 
to peak at midrapidity, suggesting a significant change 
in the longitudinal phase space distribution between the 
SPS and RHIC. Measurement of the rapidity distribution 
of identified particles is needed to establish the boost in- 
variance at RHIC. 

The h~ density at midrapidity for p± > lOOMeV/c 
is dN/dri\ri=Q = 261 ± l(stat) ± 17(syst). Extrapolation 
to pj^ = yields dN/dr)\r,=o = 280 ± l(stat) ± 20(syst). 
Assuming an average of 172 participant pairs per central 
Au-|-Au collision, this corresponds to 1.63 ± 0.12 ft," per 
participant nucleon pair per unit pseudorapidity, a 38% 
increase over the yield in pp collisions extrapolated to the 
same energy (we neglect isospin correction factors of 
order 1-3%) and a 52% increase over Pb+Pb collisions 
at ^= 17 GeV §. 

For the total charged multiplicity density in Au+Au 
interactions at ^/s^= 130 GeV, the PHOBOS collab- 
oration has reported dNch/dr]\^^^^i = 555 ± 12(stat) ± 
35(syst) (6% most central collisions) while the 

PHENIX collaboration has reported dN ch I dr]\n=o = 
622± l(stat) ±41(syst) (5% most central) To com- 
pare to these results, positive charged particles were anal- 
ysed in the framework described above. For the 5% most 
central collisions, STAR measures a total charged multi- 
plicity density of dN^h/ drj\,-i^o = 567± l(stat) ±38(syst). 

In conclusion, we find that particle production per par- 
ticipant in central Au-|-Au collisions at ^Js^^ =130 GeV 
increases by 38% relative to pp and 52% compared to nu- 
clear collisions at ^/s^= 17 GeV. The p± distribution 
is harder than that of the pp reference system for the pj_ 
region up to 2 GeV/c. Scaling of produced particle yield 
with number of participants shows a strong dependence 
on pj^, with Wounded Nucleon scaling achieved only at 
the lowest measured p^. The h~ pseudorapidity distri- 
bution is almost constant within \ri\ < 1. 
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FIG. 1. Normalized multiplicity distribution of h~ with 
p± > 100 MeV/c in Au+Au collisions at ^/s^ = 130 GeV. 
The shaded area is 5% most central collisions, selected by 
ZDC coincidence. The solid curve is the prediction from the 
HIJING model. 
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FIG. 2. Upper panel: h~ px-spectra for the 5% most cen- 
tral Au-j-Au collisions at midrapidity (jjyj < 0.1) for several 
systems. The curves are power-law fits to the data. Lower 
panel: ratio of STAR and scaled UAl px-distributions (see 
text). 
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FIG. 3. h~ pseudorapidity distribution from 5% most cen- 
tral Au+Au collisions for p± > 100 MeV/c (filled circles) and 
all p± (open squares). 
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